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ABSTRACT 

We conduct numerical experiments by evolving gaseous/stellar disks embedded in live dark matter 
halos aiming at quantifying the effect of gas spatial resolution and gas content on the bar evolution. 
Three model sequences have been constructed using different resolution, and gas fraction has been 
varied along each sequence within the range of /g= 0% — 50%, but keeping the disk and halo properties 
unchanged. We find that the spatial resolution becomes important with an increase in the gas content. 
For the higher resolution model sequences, we observe a bimodal behavior in the bar evolution with 
respect to the gas fraction, especially during the secular phase of this evolution. The switch from the 
gas-poor to gas-rich behavior is abrupt and depends on the resolution used, being reasonably confined 
to /g~ 5% — 12%. The diverging evolution has been observed in nearly all basic parameters charac- 
terizing bars, such as the bar strength, central mass concentration, bar vertical buckling amplitude, 
bar size, etc. We find that the presence of the gas component severely limits the bar growth and 
affects its pattern speed evolution. Gas-poor models display rapidly decelerating bars, while gas-rich 
models exhibit bars with constant or even slowly accelerating tumbling. We also find that the gas-rich 
models have bar corotation (CR) radii within the disk at all times, in contrast with gas-poor and 
purely stellar disks. In addition, the CR-to-bar size ratio is less than 2 for gas rich-models. Next, we 
have confirmed that the disk angular momentum within the CR remains unchanged in the gas-poor 
models, as long as the CR stays within the disk, but experiences a sharp drop before leveling off in the 
gas-rich models. Finally, we discuss a number of observed correlations between various parameters of 
simulated bars, such as between the bar sizes and the gas fractions, between the bar strength and the 
buckling amplitude, the bar strength and its size, etc. 

Subject headings: stars: kinematics and dynamics — galaxies: evolution - galaxies: halos - galaxies: 
kinematics and dynamics - galaxies: spiral — galaxies: structure 



1. INTRODUCTION 

Many issues related to the formation and evolution of 
galactic bars remain unsettled (e.g., Villa- Vargas, Shlos- 
man & Heller 2009, hereafter Paper I). In Paper I, we 
have revisited the properties of barred disks embedded in 
the dark matter (DM) halos, with an emphasis on the an- 
gular momentum redistribution in collisionless systems. 
Here we attempt to understand some effects of the gas 
component on the bar evolution. In a follow-up work, we 
investigate the effect of varying the halo parameters in 
the gaseous/stellar disk models. 

Due to its viscous and dissipative nature, the gas 
can infiuence the stellar component in the disk, well 
beyond its observed mass fraction (e.g., Shlosman & 
Noguchi 1993; Friedli & Benz 1993; Heller & Shlosman 
1994; Berentzen et al. 1998; Bournaud & Combes 2002; 
Berentzen et al. 2007; Curir et al. 2007). The gas dis- 
sipation, especially in the bar's presence, triggers the 
mass redistribution in the disk. A substantial fraction 
of barred orbits are self-intersecting, which while irrel- 
evant for stars has a strong effect on the gas. The gas 
populating such orbits is short-lived, resulting in shocks 
and loss of rotational support. Different viscosities in 
gas and stellar 'fiuid' lead to a delayed response to the 
gravitational torques and to an exchange of angular mo- 
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mentum between these components. The loss of angular 
momentum in the gas induces a flow down to the in- 
ner disk, forming a central mass concentration (CMC), 
and possibly forming and fueling the central supermas- 
sive black hole (SBH) via nested bars mechanism (e.g., 
Shlosman et al. 1989, 1990; Pfenniger & Norman 1990; 
Friedli & Martinet 1993; Knapen et al. 1995; Begelman 
& Shlosman 2009; Hopkins & Quataert 2010), and addi- 
tional processes closer to the SBH (e.g., Shlosman 1999; 
Hopkins et al. 2010). The gas is capable of modifying 
the disk orbital structure (e.g., Berentzen et al. 1998). 
Moreover, the cold gas is inherently clumpy which effec- 
tively heats up the stellar 'fiuid' and populates the disk 
with an increased fraction of chaotic orbits, resulting in 
a decay of the bar strength (e.g., Shlosman & Noguchi 
1993; Berentzen et al. 2007). The vertical buckling insta- 
bility in the bar is progressively damped by the two-fiuid 
gas-star interaction (e.g., Berentzen et al. 2007). 

Formation and evolution of stellar bars is intricately 
related to the redistribution of angular momentum, J, in 
the disk-halo system. This means the presence of sources 
and sinks of angular momentum. Sellwood (1981) has 
shown that the former reside in the disk and the latter 
in the DM halo (see also Debattista & Sellwood 1998; 
Athanassoula 2002). Interactions of sources and sinks of 
J are dominated by resonances (Lynden-Bell & Kalnajs 
1972; Athanassoula 2002; Martinez- Valpuesta, Shlosman 
& Heller 2006; Paper I). Details of these resonance inter- 
actions are not fully understood and quantified. This is 
especially true in the presence of the gas component, as 
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the above cited works are based on a purely collisionless 
modeling. 

In this paper we analyze the effect of the gaseous com- 
ponent on the evolution of galactic bars and take a nu- 
merical approach. In particular, we vary the gas content 
in the disk and advance models with various numerical 
resolution, but keep the stellar disk and DM halo para- 
maters unchanged, using the Standard Model of Paper I. 
We are especially interested in how the gas presence and 
the associated numerical resolution affect the basic pa- 
rameters of a stellar bar, as well as the angular momen- 
tum transfer process in a barred disk and in a disk-halo 
system. For this purpose, we perform a detailed com- 
parison of models with gas to the standard collisionless 
model published in Paper I. 

2. NUMERICS AND MODELING 

We use the hybrid TV-body and Smooth Particle Hy- 
drodynamics (SPH) FTM-4.4 code (e.g.. Heller & Shlos- 
man 1994; Heller, Shlosman & Athanassoula 2007; 
Romano-Diaz et al. 2009) to evolve the stellar and 
gaseous disks embedded in the DM halos. The gravi- 
tational forces are calculated using the FalcON routine 
(Dehnen 2002) which scales as 0{N). The adopted units 
are the same as in Paper I: the units of mass and dis- 
tance are taken as 10^^ Mq and 10 kpc respectively. 
This makes the unit of time equal to 4.7 x 10^ yr, when 
G = 1, and the velocity unit 208 km s^^. The gravita- 
tional softening is egrav = 0.016 for stars and DM par- 
ticles. For the gas particles we use a dynamical soften- 
ing. The gravitational softening is set to the smoothing 
length unless the smoothing length falls below the fixed 
limiting value edyn- The models consist of a stellar disk 
with A^, = 2 X 10^, a gas disk with A^gas = 4 x 10"* and 
DM halo with TYdm = 10^ collisionless particles. Models 
were evolved for about a Hubble time. At = 270 in the 
adopted units, which translates to 12.7 Gyr. 

During the model evolution, we have routinely ob- 
served the formation of very compact accumulations of 
gas particles in the central disk region. In the mod- 
els with a substantial gas component we have observed 
the formation of secondary gaseous bars which decoupled 
from the large-scale bars and contracted subsequently to 
spatial scales where insufficient resolution resulted in flat- 
tened 'blobs' a few softening lengths in radius, in agree- 
ment with simulations of Englmaier & Shlosman (2004) . 
Besides the gradual capture of additional gas particles, 
the morphological structure of the blob has evolved very 
little and was beyond the resolution of our models. For 
the sake of shortening the computational time, we have 
replaced the gas particles trapped in the center by stel- 
lar ones. This operation was repeated whenever the cen- 
tral gas accumulation was substantially slowing down the 
overall evolution. We kept the softening of the individ- 
ual particles when they were converted from one type to 
the other. We have run a large number of tests to verify 
that this action did not affect the evolution of the bar. 
This was achieved by running parallel models with and 
without the gas particle replacement. 

2.1. Initial conditions and model parameters 

The initial conditions of the stellar and DM particles 
were created with the procedures described in Paper I, 
using the density profiles from Hernquist (1993). The 



mass volume density distribution in the disk is given in 
cylindrical coordinates by 
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where Md is the disk mass, /i is a radial scale length and 
zq is a vertical scaleheight. The density of the spherical 
halo is given by 
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where is the mass of the halo, r^ is a Gaussian cutoff 
radius and 7 is the core radius, a is the normalization 
constant defined by 
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with q — 7/rc. The particle velocities, dispersion ve- 
locities and asymmetric drift corrections were calculated 
using moments of the collisionless Boltzmann equation. 
Since models thus constructed are not in exact virial 
equilibrium, the halo component was relaxed for t ^ AO 
in the frozen disk potential. 

Because we are interested to quantify the effect of the 
gas fraction and gas spatial resolution on the bar evolu- 
tion, we use the pure stellar model SD from Paper I as 
our benchmark model (Table 1). A fixed fraction /g of 
stellar disk particles at i = were converted to identical 
mass gas particles and re-balanced using the central at- 
traction forces from the total mass distribution. The gas 
is considered to be isothermal with Tgas = 10"* K, and 
initially moves on circular orbits. 

We have created a set of models covering a two dimen- 
sional parameter space: by varying the gas mass frac- 
tion, /g, in the disk, and by changing the limiting value 
of the gravitational softening, edyn, in the gas. The sum 
of stellar -|-gas mass was kept constant in the models. We 
used the values /g = 0%, 2%, 4%, 8%, 15%, 30% and 
50%, and egrav= 0.016, 0.05 and 0.1. The softening used 
for the stellar and DM components is fixed at 0.016. Ta- 
ble [5] shows the combination of /g and egrav values used 
in each model. 



TABLE 1 

Parameters of the Standard Model 



HALO 
Parameter Value 


DISK 
Parameter Value 


A^DM 


10" 


A. 


2 X 10^ 


Mh 


3.15 


A/d 


0.63 


rt 


8.55 


Rt 


1.71 


7 


0.1425 


h 


0.285 


rc 


2.85 


za 


0.057 






Q, 


1.5 



Note. — Q, is the Toomre parameter for the stellar component 
fixed at R = 2.4/i, where h is the thickness of the disk; rt and Rt 
are numerical truncation radii in the halo and the disk. All values 
are given in dimensionless units, §2. 



3. RESULTS 

Disks with a high content of gas tend to form dense 
clumps of gas which interact with the stellar component 
and raise the velocity dispersion in the disk. If this rise in 
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TABLE 2 

Gas fractions and limiting dynamical softening in the 
model sequences 



Model 




^grav 


bU 




U 




SD. 


.G2S1 




U.UlD 


SD. 


.G4S1 


4 


0.016 


SD. 


.G8S1 


8 


0.016 


SD. 


.G15S1 


10 


U.UlO 


SD. 


.G30S1 


Qn 
oU 


U.UID 


SD. 


.G50S1 


OU 


U.UID 


SD. 


.G2S2 


Z 


U.UOU 


SD. 


.G4S2 


4 


0.050 


SD. 


.G8S2 


8 


0.050 


SD. 


.G15S2 


15 


0.050 


SD. 


.G30S2 


30 


0.050 


SD. 


.G50S2 


50 


0.050 


SD. 


.G2S3 


2 


0.10 


SD. 


.G4S3 


4 


0.10 


SD. 


.G8S3 


8 


0.10 


SD. 


.G15S3 


15 


0.10 


SD. 


.G30S3 


30 


0.10 


SD. 


.G50S3 


50 


0.10 



Note. — All values are given in dimensionless units, [[2] 
Columns (1) model sequences; (2) gas fractions in %; (3) gravi- 
tational softening in the gas. 



the disk 'temperature' is substantial, the bar instability 
may be lessened or completely suppressed (Shlosman & 
Noguchi, 1993). We have encountered this problem when 
evolving models with /g=100% which prevented us from 
completing these runs. 

3.1. Bar strength 

The stellar bar strength has been quantified using the 
Fourier amplitude of the m — 2 mode normalized 
by the m = mode (Paper I). It is obtained by in- 
tegration over restricted cylindrical volumes where the 
bar is a dominant morphological feature, namely, over 
R = 0.1— i?b range, where i?b is the bar size defined 
in Section 3.3 (see also Paper I). In Figure [U we plot 
as a function of time t for various gas fractions, /g, 
and gravitational softening in the gas, egrav The initial 
stages in the evolution of are very similar in all mod- 
els: an initial stage of an accelerated growth, and a peak 
followed by a sudden drop. The duration of this dynam- 
ical stage varies from model to model, but is completed 
by t - 100. 

The first peak in is followed by the vertical buck- 
ling instability in the bar leading to an abrupt weakening 
in the bar but not a complete dissolution (e.g., Martinez- 
Valpuesta & Shlosman 2004). We observe that the bar 
weakening is quite independent of fg and Cgrav, as ex- 
pected. For gas-poor models, the first peak of A2h is 
independent of egrav and fg. For gas- rich models, with 
/g^ 15%, the peak is lowered gradually for Cgrav^ 0.05, 
up to a factor of 0.4. Models with egrav= 0.1 appear not 
to be affected at all by this trend. 

The post-buckling evolution of the bars is much more 
affected by fg and Cgrav, and shows a bimodal behavior. 



Namely, in gas-poor models, the bar resumes its growth 
but at a more gradual pace as compared to the dynamical 
growth. In gas-rich models, the bar strength declines 
over the simulation (i.e., Hubble) time. This bimodal 
behavior was noted by Berentzen et al. (2007), but the 
current models show it more explicitly. 

We shall refer to the first bar evolution phase, includ- 
ing the buckling, as the dynamical phase, and to the 
subsequent evolution as the secular phase. The secular 
stage allows the separation of the secularly growing mod- 
els from the secularly declining ones. In no models have 
the bars disappeared completely — at least a substantial 
oval distortion remained. 

The borderline fg between these two trends, growth 
and decline, depends on Egrav For egrav= 0.016, it lies in 
the range of fg^ 5% — 7%. Larger Cgrav moves the border- 
line fg toward more gas-rich models. Models SD_G8S2 
and SD.G15S3 lie close to this borderline, /g- 8%- 10% 
for egrav= 0.05 and /g- 10% - 12% for Cgrav^ 0.1, and 
show a very singular evolution — after the sudden drop 
^2b~const. for about At ~ 100, and a second period of 
bar growth begins. This is a surprisingly mixed behavior 
showing the prolonged constancy in A2h of the gas-rich 
models and the secular bar growth of the gas-poor ones. 

Model SD_G50S2 has a peculiar evolution in the dy- 
namical stage that deserves special attention. Even 
though this is a gas-rich model, the bar resumes is sec- 
ular growth after the buckling. The rate of growth is 
much more moderate than that in the other models but 
is clearly noticeable. 

3.2. Pattern speed 

The evolution of the pattern speed fib is shown in 
Fig. [2] It remains about constant during the dynami- 
cal phase of the bar evolution, slightly rising in the gas 
rich models. The bimodality of gas-poor and rich models 
shows up in the secular evolution where flh drops in the 
former and stays constant or rises in the latter disks. In 
the secular phase, fib nearly always anticorrelates with 
A2h A sustained drop in i^b corresponds to a rise in A2h- 
Thus, in the gas-poor models, the bar tumbling slows 
down, while in the gas-rich models it stays nearly con- 
stant. 

Fig. 2 provides a hint to what one can expect for the 
evolution of the angular momentum in the disk-halo sys- 
tem, although in no way does this figure account for all 
the angular momentum in the disk, and even in the bar 
region it represents only the tumbling of the bar. The 
gas-poor disks lose their tumbling angular momentum 
efficiently, while gas-rich ones nearly conserve it or even 
speed up their tumbling during the secular phase. 

3.3. Bar length 

The length of the bar semi-major axis, i?b, has been 
taken as the radius where the bar equatorial ellipticity 
drops by 15% off its peak (Paper I). This method has 
been tested in comparison with alternative method based 
on the last stable orbit supporting the bar (Martinez- 
Valpuesta et al. 2006). It is the most reliable when ap- 
plied after the first maximum of the bar strength. The 
ellipticity of the bar at different radii is obtained by fit- 
ting ellipses to the isodensity contours in the face-on disk. 
As in the collisionless models, the bar length exhibits 
the initial period of an accelerated growth and reaches 
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Fig. 1. — Time evolution of normalized bar amplitudes A2h- Rows correspond to models of equal dynamical softening in the gas, egrav, 
indicated on the right. Columns correspond to models with equal gas fraction, /g, indicated on the top. All the data has been smoothed 
with a lowpass Fourier filter. 
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Fig. 2. — Time evolution of the bar pattern speed Q^,. The distribution of models in rows and columns is that of Fig.[ll Model SD_G50S2 
curve has a gap around the time of the buckling. The bar weakens abruptly and its figure becomes distorted which makes it difficult to 
determine the pattern speed. Note that for the gas-rich disks, the pattern speed is a non- decreasing function of t. The gas fractions and 
spatial resolution are indicated at the top and on the right respectively. 



a maximum which always coincides (in time) with the bar instability itself. Similarly to the pure stellar model, 
maximum in (Fig. IS]). This symbolizes the period of there is a subsequent drop in i?b as a result of the verti- 
a dynamic growth of the bar, or growth related to the cal buckling (section 3.4). During the secular phase, the 
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Fig. 3. — Time evolution of the bar semimajor-axis length f?b- The distribution of models in rows and columns is as in Fig. [T] The bar 
sizes include the ansae (see text). Note that bars are not growing in the secular phase in the gas rich disks (see also SN93), and this is 
more pronounced in high resolution simulations. Gas fractions and resolution are indicated at the top and on the right respectively. 



bar length exhibits a more complicated behavior than in 
the coUisionless models. Most importantly, it grows in 
gas poor models while it stagnates or even shortens in 
the gas rich disks. 

In the late stages of secular evolution, the gas poor 
models show a sudden decrease in the bar size which has 
no immediate correspondence to ^2b evolution. Rather 
generally, the evolution of i?b does not go in tandem with 
changes in ^2b- As an example, in model SD_G4S1 after 
t ^ 220, i?b continues its growth while growth satu- 
rates. This behavior is caused by the formation, growth 
and detachment of the ansae^, as discussed in Paper I. 
In the following, we shall show that i?b correlates with 
other properties of the bar-disk system. 

The bar corotation (CR) radius i?cr has been computed 
using linear approximation. i?cr grows or shrinks as a 
consequence of the variation of the bar tumbling speed, 
i7b- It grows substantially in the SD and gas poor mod- 
els, while it stays constant or even drops slightly in the 
gas rich models. Its initial value among all models is the 
same. As in Paper I, we follow the ratios of i?cr to that of 
the disk (Fig. SJ and bar i?b sizes. Both and i?b 
have been defined in Paper I. Clearly, when R„/Rd^ 1, 
its growth experiences a temporary slowdown and then 
resumes, rising to larger values. This effect is the re- 
sult of the combination of a temporary slowdown of the 
outwards radial movement of i?cr and a temporary ex- 
pansion of the disk radial border which results from the 

^ The ansae (i.e., handles) are typically found in early-type 
disk galaxies (Fig. 11 in Martinez- Valpuesta et al. 2006; see 
also Martinez- Valpuesta, Knapen & Buta 2007), e.g., NGC 4262, 
NGC 2859, and NGC 2950 (Sandage 1961), NGC 4151 (Mundell 
& Shone 1999), and ESQ 509-98 (Buta ot al. 1998). 



drag on some of the disk material trapped at the CR 
resonance. As in the pure stellar models, the growth of 
the bar strength is sensitive to the moment at which i?cr 
crosses the edge of the disk, seen as a temporary period 
of slower growth of A2h- Note, that Rcr/Rd is always less 
than unity in the gas rich models and it seems that for 
a fixed /g, the CR-to-disk size ratio is damped stronger 
for small egrav, i-e., for higher numerical resolution. 

The ratio Rcr/Rh shows a pronounced dip during the 
initial growth of the bar. This is followed by an equally 
abrupt rise due to the bar shortening that results from 
the vertical buckling instability. After these variations, 
both i?b and i?cr grow in tandem, and Rcr/Rb^ 1-5 — 
2 ± 0.15 in most of the models. Exceptions are models 
SD_G4S2 which has a rise and drop resulting from a tem- 
porary stall of Rb, and models SD_G30S3 and SD_G50S3 
where i?cr~ const, while the bar shrinks. Finally, a pro- 
nounced rise above Rc^/Rb— 2 occurs in models due to 
detachment of the ansae (Paper I). This rise is absent in 
models, in which no ansae are formed or which form too 
late in the run and do not have time to detach from the 
bar, as in SD_G15S3. We find that the ansae do not form 
in models with more than ^ 10% of the gas. We note 
that the bar size evolution given in Fig. 3, includes the 
ansae. However, in Section 4, we find it advantageous to 
exclude the ansae when discussing some correlations. 

3.4. Vertical buckling 

The vertical buckling observed in many numerical bars 
is an event that in some cases can reshape the phase- 
space density of the disk. As we found in Paper I, when 
the disk is heated vigorously by the buckling instability, 
the secular growth of the bar can be seriously dimin- 
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Fig. 4. — Time evolution of the ratio Rcr/Rd- The horizontal line has been added to help identify the time i?cr crosses the disk border. 
The distribution of models in rows and columns is as in Fig. [l] Note that for the gas-rich disks the CR is always inside the disk, especially 
for higher resolution disks. Gas fractions and resolution arc indicated at the top and on the right respectively. 



ished or completely halted. We have measured the ver- 
tical asymmetry of the stellar disk with the index ^i,z, 
computed as described in Paper I. The evolution of v4i ^ 
is shown in Fig. [5j Most of the models have at least one 
clear peak between t = 50 and 100, and this peak coin- 
cides in time with the drop in A2h- Exceptions are mod- 
els SD_G15S1, and SD_G50S2, which have no clear peak 
above the noise level. Some models have secondary and 
even higher order peaks or prolonged periods where Ai z 
is clearly above the noise level, as found by Martinez- 
Valpuesta et al. (2006). 

Clearly, the height of the first peak varies with the gas 
content of the disk. In the sequences with egrav= 0.016 
and 0.05, the peak becomes gradually lower in models 
with higher /g. This result is in good agreement with 
Berentzen et al. (2007). The big surprise is in the se- 
quence with egrav= 0.1, where the value of the peak drops 
as /g increases from 0% to 8%, and then rises in mod- 
els with /g= 30% and 50%. This is discussed further in 
Section 111 

The signature of the vertical buckling can be observed 
in the behavior of J7b- As the bar strengthens, it brakes 
against the outer disk and against the DM halo, at later 
time and 57b declines. At the time of buckling, the bar 
weakens abruptly and flh experiences a break — its slope 
changes and becomes much smaller. 

3.5. Central mass concentration 

One of the differences that arise between pure stellar 
models of galactic bars and models including a gas com- 
ponent is the formation of a central mass concentration. 
The central densities, both stellar and DM, have notably 
increased even in pure stellar models as a response to the 



asymmetric potential of the bar and the vertical buckling 
(Paper I; Dubinski, Berentzen & Shlosman 2009). How- 
ever, this effect is by far more intense in models with 
gas. The growth of central mass concentration (CMC) 
is boosted by the bar which channels the gas toward the 
central kpc. This leads to a CMC composed mainly of 
gas, and, in smaller proportion, of stellar and DM par- 
ticles. In all our models, the CMC has resided inside a 
radius smaller than 0.1. 

To follow the growth of the CMC, we have measured 
the mass Mcmc contained within a sphere of radius 0.1 
placed at the center of mass of the disk. Contributions 
from all three components have been included (Fig. 6). 

Model SD_G2S1 serves as a benchmark of evolution in 
the gas poor models. With the formation of the stellar 
bar, a massive accumulation of gas appears at the center, 
with an elongated shape of ^ 0.1 x 0.04, in the disk plane. 
This CMC is approximately aligned with the major axis 
of the (gas) bar. It rapidly captures the gas fed by the 
bar and simultaneously contracts into a small, more ax- 
isymmetric and somewhat flattened circular 'blob' with 
AR ^ 0.04 and Az ^ 0.01. These characteristic sizes are 
determined largely by Cgrav The accelerated Mqmc rise 
seen in Fig. 6 encompasses the exponential bar growth 
phase, up to the time of buckling (or the time when 
drops in models with no buckling). Typically, about 50% 
of the gas mass is captured by the CMC during this stage. 

During the secular phase of the bar evolution, the CMC 
captures the gas at a much slower pace. This is mostly 
the gas which remains outside i?b and even _/?„• at the 
onset of the secular phase. The slowing down of the CMC 
growth results from a combination of two factors. First, 
the availability of the disk gas has been severely reduced 
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Fig. 5. — Evolution of the bar vertical asymmetry measured by the Fourier coefficient |j4i_z| of the m = I mode in the rz-plane corotating 
with the bar major axis. No filtering was applied to this data. The distribution of models in rows and columns is as in Fig. [l] Note that 
the first buckling amplitude decreases with /g, but the sequence egrav= 0.1 exhibits a somewhat more complicated behavior. In addition, 
Ai.zl correlates with egrav Gas fractions and resolution are indicated at the top and on the right respectively. 
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Fig. 6. — Time evolution of the central mass concentration, CMC: the total mass contained inside a sphere of radius 0.1 at the center of 
mass of the stellar disk. The stellar, gaseous and DM components are included in the mass count. All the data has been smoothed with a 
high frequency Fourier filter. Gas fractions and resolution are indicated at the top and on the right respectively. 



by the violent initial inflow, and second, the bar strength in the disk after the initial inflow, although they show 
has been diminished by the buckling. As a consequence, a healthy secular bar growth. Whereas gas-rich models 
the gas-poor models have been left with very little gas have gas left, but their bars are weak and do not grow. 
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Fig. 7. — Time evolution of disk total angular momentum, J^. Gas fractions and resolution are indicated at the top and on the right 
respectively. 



SO no fresh gas crosses Ra which stagnates as well. 

Three models exhibit a somewhat different evolution 
for AfcMC- Instead of a short period of a very acceler- 
ated growth, the CMC grows at a somewhat slower rate 
but over an extended period of time. These are mod- 
els SD_G15S3, SD_G30S3 and SD_G50S3, which are all 
gas-rich and have the lowest resolution. In these models, 
the gas content of the bar has formed without a visibly 
prominent CMC. Instead of a rapid influx to the center, 
the gas in the bar contracts rather slowly, gradually in- 
creasing its density at the center. The process saturates 
when the gas is completely contained within the central 
Rr^ 0.1. 

Models SD_G8S2 and SD_G15S3 show a rate of Mcmc 
growth that increases with time at the end of the run. 
This behavior is related to the late period of accelerated 
bar growth. 

3.6. Angular Momentum 

Angular momentum evolution in the disk has been fol- 
lowed within a number of characteristic radii, namely, the 
CR ( Jd,in) and the disk radius, (which contains 98% of 
the disk mass by definition). The latter exhibits a clear 
evolutionary sequence as a function of /g (e.g.. Fig. 7). 
Both the amount of J lost by the disk over dynamical 
and secular evolution decreases monotonically from pure 
stellar disks to progressively more gas-rich ones. In a 
way this is a reflection of weaker bars along /g. Even 
more interesting is the evolution of Jd.in (Fig- 8). After 
the initial adjustment, the disk within the CR is los- 
ing its angular momentum. The loss is not dramatic, 
^ 10% — 20%, and somewhat increases with /g. This de- 
cline in Jd,in continues until the vertical buckling sets in. 
For low /g disks, the buckling is associated with a nearly 



sudden increase in Jd,in which restores the pre-buckling 

value followed by a subsequent gradual decline. Gas-rich 
disk show no sudden increase in Jd,in but rather a slow 
decline. 

Jd,in stays quite flat until t ~ 150, a long time af- 
ter the buckling, as can bee seen from Fig. 8 (with the 
exception of a small initial decline and increase, as men- 
tioned above). Interestingly, this time corresponds to 
Rci/Rd< 1 in the gas-poor models (Fig. 4). Crossing 
Rcr/Rd= 1 border affects the bar growth (Fig. 3), which 
saturates immediately. It also is reflected in the evolution 
of Rci/Rb which increases abruptly thereafter. Clearly, 
the bar growth ceases at some point when the Rcr ex- 
ceeds _Rd and the bar cannot capture additional orbits 
and be fed by the angular momentum from these orbits. 
This situation is similar to that analyzed in Paper I (Sec- 
tion 3.6), where the angular momentum of the disk inside 
i?cr stays about constant as long as R^ remains within 
the disk. Our present models can be directly compared 
to the Standard model of Paper I. We return to this point 
in Section 4. 

The loss of the angular momentum, Ja.out) with /g 
by the outer disk appears to be much more severe than 
Jd.in- But one should remember that the outer disk, 
i.e., outside the CR, accounts for less mass. Before the 
buckling, Jd,out behaves similarly to Jd,in, but differently 
at the later times. The trend here is that Jd,out declines 
steeply with time for gas poor disks after the buckling 
period. It stays nearly constant with time for the gas 
rich models. 

Models with different Cgrav show a behavior consistent 
with our understanding. There is a slight increase in the 
J transfer during dynamical and secular stages of the 
bar evolution. Much less difference is observed in the 
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Fig. 8. — Time evolution of disk angular momentum inside the bar corotation, Jd.in- Note that Jj^inRiconst. in time as long as Rcr 
remains inside the disk, as can be seen from Fig. 4. Gas fractions and resolution are indicated at the top and on the right respectively. 



evolution of Jd,in and Jd,out in this case. 

4. DISCUSSION: TESTING NEW CORRELATIONS 

Wc have studied the effects of gas fraction, /g, and spa- 
tial resolution, Cgrav, on some aspects of bar evolution in 
galactic disks embedded in DM halos. Specifically, we 
aimed at understanding their effect on the basic param- 
eters characterizing stellar bars, e.g., bar size, corotation 
radius, strength, etc. We have also followed the angular 
momentum redistribution in the disk-halo systems as a 
function of /g and Cgrav The SD pure stellar model (Pa- 
per I) has been used as a template — the disk and halo 
parameters in current models have been fixed at those of 
the SD. As before, the bar evolution has been separated 
into two main phases — the dynamical phase, which rep- 
resents the bar instability itself de facto, and the secular 
phase. These phases arc separated by the vertical buck- 
ling instability in the bar. We investigate a number of 
new correlations which are corollaries to evolution dis- 
cussed in section 3. In the follow-up work (Villa- Vargas 
et al., in preparation), we analyze the angular momen- 
tum transfer and bar evolution as a function of disk gas 
fraction and by varying the basic parameters of the DM 
halo. 

Overall, there appears to be a substantial difference 
in the evolution of the gas-poor and rich models, as has 

been shown in the previous section. This is discussed 
below in more detail. The boundary between the gas- 
poor and gas-rich models depends on the spatial resolu- 
tion for the gas component — it lies around 5%-7% for 
egrav=0.016, and shifts to ~ 10% - 12% for egrav=0.1. 
The loss of spatial resolution, in a way, changes the na- 
ture of the gas component — it becomes less dissipative. 
We define the bar strengthening during its dynamical 



phase as AA2b =A2\y{t = tpeak) -Mh{t = 0), where 
ipeak is the time when the bar has reached its maxi- 
mal amplitude prior to the onset of buckling; the second 
term A2h{t = 0)= 0. Similarly, the bar growth during 
the secular phase is defined as A^2b =^2b(^ = 270) - 
^2b(i = train); wlicrc t-cain IS thc time of tlic minimal bar 
amplitude immediately following the buckling. 

In the dynamical stage of the bar evolution, the largest 
differences have been observed in the gas-rich models, 
where AA2b drops dramatically with /g for egrav= 0.016 
and 0.05, while it exhibits no dependence on /g for 
egrav= 0.1 (Figs. 1, 9a). In the gas-poor models, the 
maximal bar amplitude achieved in the dynamical phase 
is nearly independent of the gas fraction, nor does it de- 
pends on the spatial resolution used in the modeling of 
the gas component (Figs. 1, 9a). Fig. 9a supports the 
Berentzen et al. (1998, 2007) analysis that the increas- 
ing presence of the gas component makes the bar insta- 
bility milder. In the secular stage, this sharp decrease 
in AA2b with /g persists, including the lowest resolution 
sequence as well. The bar growth is severely restrained 
in the gas-rich models in this phase. It either had van- 
ished or even becomes negative. The underlying physical 
process has been quantified by Berentzen et al. (1998): 
the gas affects the orbital structure of the coUisionless 
components, which in turn modifies the bar properties. 
However detailed understanding of how this influences 
particular parameters of the bars requires more work. 

We have attempted to clarify this dependence of the 
bar strength on /g and egrav As a working hypothe- 
sis, we have tested whether the final mass of the CMC, 
McMC, in each model can serve as an underlying hid- 
den parameter replacing /g and Cgrav The dependency 
of the final bar strength on the CMC mass in the pure 
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Fig. 9. — Bar strengthening during (a) dynamical and (b) secular phases of its evolution as a function of the gas fraction, /g, for different 
spatial resolution in the gas, egrav, namely egrav= 0.016 (crosses, solid lane), 0.05 (x, short dashed lane), and 0.1 (circles, long dashed 
lane). The filled square represents the SD model (Paper I) with no gas. The model SD_G8S1 has been omitted (see the text). 



stellar disks has been analyzed by Athanassoula, Lam- 
bert & Dehnen (2005). Extensive study of the bar evo- 
lution in the presence of an analytical CMC has been 
performed by Shen & Sellwood (2004). Bornaud et al. 
(2005) have attributed the bar weakening to both the 
CMC and the angular momentum transfer from the gas 
to the stellar bar, using analytical CMC and DM halo. 
However, the latter process has been found not important 
by Berentzen et al. (2007) following a detailed analysis of 
the angular momentum transfer in live potentials of the 
CMC and DM halo. Moreover, Berentzen et al. (1998, 
2007) have used self-consistently growing CMCs in the 
gaseous/stellar disks to demonstrate that the secular bar 
growth is strongly affected. The set of numerical models 
analyzed here is in fact the most controlled experiment 
performed so far to test the influence of gas fraction and 
its resolution on the bar evohition. 

Fig. 10a shows the dependence of the bar strengthen- 
ing after the first buckling, AA2h on final Mcmc, with 
all other parameters characterizing the stellar disk and 
the DM halo being fixed. One model, SD_G8S1, has been 
omitted, as its CMC secular evolution is unusually flat 
and falls out of the sequence with egrav= 0.016. If Mcmc 
is the single defining variable which controls AA2h, all 
three egrav-scquences are expected to merge into a single 
sequence in Fig. 10a. What is actually observed is that 
the three sequences exhibit a very similar behavior — 
each curve is flat for the gas-poor models, experiences an 
abrupt drop, and flattens out. However, the curves do 
not coincide completely as would be expected if Mcmc 
were the unique underlying parameter. In fact, in some 
cases the points corresponding to the same AfcMC but 
different Cgrav have a substantial vertical dispersion in 
AA2b- Hence, the issue remains inconclusive. In con- 
trast, we show the final Mcmc value normalized by the 
disk mass at t = as a function of /g (Fig. 10b). This 
fractional Mcmc value is indeed unique for the three se- 
quences, as all three curves have nearly merged. 



We now turn to a plausible correlation between the 
buckling amplitude and the bar strengthening. The ver- 
tical buckling (e.g., Toomre 1966; Combes et al. 1990; 
Raha et al. 1991) is a recurrent instability (Martinez- 
Valpuesta et al. 2006) that was recently analyzed by 
Berentzen et al. (2007) in the presence of gas (see also 
Berentzen et al. 1998). The gas component, it has been 
concluded, leads to a milder instability. Here we have 
attempted to relate the buckling amplitude, ^i,z, to the 
change in the bar amplitude, AA2h, in the dynamical 
and secular phases of evolution (Fig. 11). In both phases 
we observe a clear correlation between Ai_z and Aj42b. 
Stronger bar instability leads to a stronger buckling, and 
indeed, increasing /g makes the bar and buckling in- 
stabilities milder (Figs. 11a; see also Fig. 9a). On the 
other hand, stronger buckling goes in tandem with the 
bar secular growth (Fig. lib). This trend shows satura- 
tion for the strongest bars. Higher resolution sequences 
with egrav= 0.016 and 0.5 behave in a very similar fash- 
ion, while the lowest resolution sequence stands out of 
this correlation. We have also checked the value of the 
drop in A2b during the buckling and find a clear match 
between the higher resolution sequences, while the lower 
resolution models behave differently. 

The physical extent of the bar depends on its ability to 
capture additional orbits. While in principle this capture 
can proceed at all radii, the fertile region lies between the 
bar's end and its CR radius, where various families of or- 
bits can be easily destabilized. Therefere, the bar growth 
due to the orbit capture should go in tandem with the an- 
gular momentum influx because the near-CR orbits will 
have a larger momentum-to-energy, J/E, ratio than the 
bar orbits. In Paper 1, we have shown, and this is con- 
firmed here, that the influx of angular momentum across 
the CR and into the bar is able to maintain Jd,in~const. 
in time, as long the CR radius lies within the disk, but 
in the presence of the gas this statement is limited to the 
gas-poor models only. This happens despite that the CR 
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Fig. 10. — (a). The secular strengthening of the bar, AA2b, after the first buckling and until the end of the simulations, as a function 
of the final central mass concentration CMC within R = 0.1, McmCi for different spatial resolution in the gas, abbreviated as in Fig. 9. 
The filled square represents the SD model (Paper I) with no gas. The model SD_G8S1 has been omitted (see the text), (b). Final Mcmc 
normalized by the total mass of the disk at t = as a function of /g. 




Fig. 11. — First vertical buckling amplitude of the stellar bar, Ai.z, as a function of the bar strengthening during dynamical and secular 
phases of its evolution, AA2h, for different spatial resolution in the gas, abbreviated as in Fig. 9. The filled square represents the SD model 
with no gas. 



remains within the disk at all times and for all gas-rich 
models (Fig. 4). In other words, in gas-rich models the 
J influx across the CR cannot compensate for its loss 
within the CR. 

The amplitude A2h is expected to be related to the 
ability of the bar to capture additional orbits in the bar 
CR region, leading to its geometrical growth. In both 
the dynamical and secular phases of the bar evolution we 
find that the final bar amplitude, ^2b, correlates with the 
final bar size^ (Fig. 12). So stronger bars appear to be 
longer as well. Taken at face value, this property of bar 

^ The only exception appears to be the low resolution sequence 
in the dynamical stage whose points cluster in the same region 
(Fig. 12a) 



evolution appears to be supported by recent observations 
in that the K-hand images of barred galaxies exhibit a 
correlation between the bar amplitude, measured by the 
m = 2 Fourier component, and its size (Elmegreen et al. 
2007). However, there are caveats. 

We find that the bar sizes are substantially smaller 
in the gas-rich models compared to the gas-poor ones 
(Figs. 3, 13). While the median bar size at the end of 
the simulation in the latter models appears to be ~ 1.2, 
it is around 0.5 in the former ones, i.e., 12 kpc vs 5 kpc, 
embedded in the identical disks and halos. This differ- 
ence arises almost entirely during the secular phase of 
the bar evolution. More precisely, by the end of the dy- 
namical stage, one can notice that the bar size slightly 
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^ 0.4 0.5 0.6 0.7 0.8 0.2 0.4 0.6 0.8 1.0 1.2 

R (peak) - no ansae R (t=270) - no ansae 



Fic;. 12. — The bar amplitude, A2b, at its pre-buckliiig peak (left) and at the end of the simulation (right) as a function of the bar size 
(radius), taken at the same time, during dynamical and secular phases of the bar evolution. Different spatial resolution in the gas has been 
abbreviated as in Fig. 9. The full square represents the SD model with no gas. We have subtracted the ansae from the bar sizes given in 
Fig. 3 (more details in the text). 
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Fig. 13. — Bimodal distribution of the final bar sizes, Rh{t = 
270), as a function of the gas fraction, /g. We have subtracted the 
ansae from the bar sizes given in Fig. 3 (more details in the text). 



anti-correlates with /g for the gas-rich models and only 
for egrav= 0.016 and 0.05 sequences. This trend is suffi- 
ciently weak and bars differ by not more than 10%-15% 
of their length. So essentially, the bar growth in the dy- 
namical phase is independent of /g. 

What is rather striking, is the abrupt change in the 
bar's growth 'habit' in the secular phase — bars grow in 
the gas-poor models and stagnate in the gas-rich ones. 
This leads to a bimodal evolution of the final bar sizes 
with respect to the gas mass fraction (Fig. 13). We note 
that this diverging evolution should be observable only 
in the long-lived bars which naturally reside in the older 
disks. The relevant evolutionary timcscale is that of a 
few Gyrs. The Milky Way disk is a good candidate, as 
it is probably about 10 Gyr old and was not affected by 



major mergers over this time period (e.g., Gilmore, Wyse 
& Norris 2002). 

The anti-correlation between the final bar sizes and 
disk gas fractions, found here, complicates the simplis- 
tic picture of bar evolution. The reason for this is 
that our gas-rich models lead to more massive CMCs 
which increase bulge-to-disk mass ratios and, therefore, 
should be associated with earlier type disks. This means 
that smaller bars in our models lie in early-type disks 
(Fig. 13), while they are expected to be found in late- 
type disks, as noted by a number of surveys (e.g., Erwin 
2005; see also Laine et al. 2002). The simplest resolution 
of this discrepancy can lie in the omission of star forma- 
tion processes and especially the feedback from stellar 
evolution and from the central supermassivc black holes 
in this work. This leads to a gross over-estimate of the 
amount of gas which reaches the central kpc and hence 
contributes to the growth of the CMC by dragging stars 
and even DM inwards. We note, however, that the pur- 
pose of this numerical exercise was to understand the 
effect of gas fraction and gas spatial resolution on the 
disk-bar-halo interactions in the system. Therefore, we 
have simplified the long list of processes known to affect 
the galaxy evolution at large. 

In summary, we find that the spatial resolution in the 
gas component becomes increasingly important for the 
bar evolution in the gas-rich disks. This is true for the 
dynamical but especially for the secular phase of evolu- 
tion. In most cases, model sequences with egrav= 0.016 
and 0.05 show a similar behavior, while differing substan- 
tially from the sequence with egra.v= 0.1. 

A bimodal behavior has been fotmd for models based 
on their gas fractions. The border line between the gas- 
poor and gas-rich systems appears to lie around 5% — 
7% for higher resolution models. It shifts to ^ 10% — 
12% for the lowest resolution models. The switch from a 
gas-poor to a gas-rich behavior appears to be sufficiently 
abrupt. It is clearly visible in all basic characteristics of 
bar evolution, such as the bar strength, the CMC mass, 
the bar buckling amplitude, the bar size, etc. The largest 
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differences in the evolution have been found in the secular 
phase. 

We find that the presence of the gas component 
severely limits the bar growth and affects its pattern 
speed evolution. While pure stellar models (Paper I) ex- 
hibit a rapid slowdown of the bar tumbling, as known for 
a long time (e.g., Debattista & Sellwood 1998; Athanas- 
soula 2003), the addition of a substantial amount of gas 
reverses this trend completely. Furthermore, the CR- 
to-disk size ratio, Rc^/Rd, was determined to be an im- 
portant dynamic discriminator between various phases 
of barred disk evolution. Here we find that the gas- 
rich models are characterized by Rcr/Rd < 1 and by 
i7b ^ const. In these models, fib can even slightly in- 
crease with time. In addition, the gas-rich models main- 
tain R„/Rh < 2 which is much more in agreement with 
observations of stellar bars being fast rotators. 

The angular momentum evolution displays the same 
degree of bimodality with the gas fraction. For low /g, 
the total disk J drops steeply and monotonically with 
time after the buckling, while it decreases weakly for the 
gas-rich models. The reason for this behavior is of course 
that the bar amplitude is substantially lower in the gas- 
rich models. Our attempt to explain this difference be- 
tween the gas-poor and rich models in terms of the more 
massive CMCs in the latter ones has been rather incon- 
clusive. We shall return to this issue in the forthcoming 



work. 

Next, we have confirmed our previous claim (Paper I) 
that the angular momentum, Jd,in, within the CR radius 
is maintained at a constant level due to the influx of 
angular momentum across the CR as it expands. We 
have extended this statement to the gas-poor models. 
For the gas-rich models, Jd,in drops abruptly to a lower 
level and stays constant thereafter. 

A number of corollaries follow from the above results. 
We find that the bar strength inversely correlates with 
the gas fraction, both in the dynamical and secular 
phases of bar evolution. The only exception seems to be 
the lowest resolution sequence in the dynamical phase 
which fails to capture this trend. We also find that the 
buckling amplitude becomes larger for stronger bars prior 
to the onset of buckling. On the other hand, the secu- 
lar growth is most prominent in bars which show a large 
buckling amplitude. Finally, we show that stronger bars 
are also the longest ones throughout both evolutionary 
phases, and that bar sizes anti-correlate with the gas 
fraction. 
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